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A remark on the local density approximation with the gradient corrections and the 

X a method 

Y.-B. XieQ B.-H. Wang, and W.-X. Wang 
Department of Modern Physics, University of Science and Technology of China, Hefei, 230026, P.R. China 

We report that the solids with narrow valence bands cannot be described by the local density 
approximation with the gradient corrections in the density functional theory as well as the X a 
method. In particular, in the case of completely filled valence bands, the work function is significantly 
underestimated by these methods for such types of solids. Also, we figured out that these deficiencies 
cannot be cured by the self-interaction-corrected-local-density-approximation method. 

PACS numbers: 71.15.Mb, 31.15.Ew 



I. INTRODUCTION 

The density functional theory 1 , 0] with the general- 
ized gradient approximations [3, [j] (GGA) and the X a 
theory [j| [f| have been quite successful in the calcula- 
tion of various properties in many systems Q, such as 
the bond length, the bond angle, various properties re- 
lated to the charge density and the vibration frequencies 
of many moleculars and solids. Moreover, the binding 
energies of many moleculars and the cohesive energies of 
many solids can be as well roughly evaluated by the GGA 
and the Xq, theory. 

Despite those successes, it is also well knownQ that the 
single particle electronic structure cannot be very well de- 
scribed by the GGA in many solids, since the GGA tends 
to underestimate the valence band width and the band 
gap. This is not difficult to understand because even 
in the electron gas model, the exchange energies for the 
states far below or far above the Fermi surface are sig- 
nificantly underestimated or overestimated, respectively. 
Nevertheless, the exchange energy for the states at the 
Fermi level is correctly given by the GGA for the electron 
gas case. Therefore, the properties related to the elec- 
tronic structure at the Fermi surface, for instance, the 
work function, has been considered correctly explained 
by the GGA for many solids. 

On the other hand, the eigenvalues of the X a single 
electron equation cannot be used to estimate the ioniza- 
tion energies of atoms and small moleculars. Fortunately, 
this problem can be cured by introducing the so called 
transition state method [5|. 

Let us consider the solids in which the valence elec- 
tronic states deviates from the plane waves significantly. 
The electronic properties of this type of solids cannot be 
described by the electron gas model. In contrast, some 
characters of the individual atom or small molecular are 
still kept in such solids. But since the single particle 
states are quite extended in crystalline solids, the tran- 
sition state method will no longer be helpful. Hence, a 
natural question is addressed: can the GGA or the X a 



method correctly describe the single particle states near 
the Fermi surface for such type of solids? We figured out 
that the answer to this question, unfortunately, is nega- 
tive from both the academic and practical point of view. 
Our paper will provide enough evidence to support this 
conclusion. 

In this paper, we study a system composed of spin 
polarized hydrogen atoms being located on a cubic lat- 
tice. The reason to choose such system is for simplicity. 
Since the spins of all electrons are parallel, the correla- 
tion energy is insignificant, and actually is zero in the 
case of large lattice constant. The exchange energy can 
be treated exactly by the Hartree Fock method. Also the 
exchange functional in the GGA and the X Q method is 
simple and the correlation functional is negligibly small 
in this system. 

In section III, we consider the large a case in which 
the overlap of the localized hydrogen electron wavefunc- 
tions on different sites is negligible. We rigorously show 
that in GGA or in X a method, the ionization energy in 
such case is given by the eigenvalue of the single particle 
equation for one hydrogen atom which is roughly[9j 0.29 
and deviates from the exact ionization energy 0.5 — 37/a 4 
considerably (Note that here the term 37/a 4 is the hole 
polarization term.). This difference, as we can see, is 
present not only in the spin polarized hydrogen solids, 
but also in many other systems, where atoms are well 
separated and the single particle wavefunctions are ex- 
tended. Notice here that the total bulk energy can be 
still correctly given by the GGA or the X a theory. 

In section IV, we study the finite a case in which the 
overlap of the localized wavefunctions on different sites 
cannot be neglected. In this case, the single particle 
wavefunctions must be extended. Therefore, similar to 
the reported results in the last paragraph, we expect that 
the work function is still significantly underestimated by 
the GGA or the X a method when a is not small. This ex- 
pectation is indeed correct up to a = 5 at which the band 
width is nearly 0.25. When a is sufficiently small, the sin- 
gle particle states can be characterized by the plane wave 
and the work function should be correctly predicted by 
the GGA or the X Q theory. 
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Finally in section V, we summarize our results. 
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II. THE MODEL 

Let us consider N hydrogen atoms located on a cubic 
lattice with the lattice constant a. For simplicity, we 
assume that the spins of all electrons are up. In this 
paper, we study the work function for the N electrons 
system, i.e., the ground state energy difference between 
the N electrons system and the N — I electrons system 
with fixed a. 

We investigate this system using the Hartree Fock 
method and the GGA or the X Q method. The corre- 
lation energy is not significant because the spins of the 
electrons are parallel. Actually, the contribution of the 
correlation energy to the work function is estimated to 
be — 37/a 4 when a is large. Hence, the Hartree Fock 
method is accurate to predict the work function. On the 
other hand, the GGA is also quite simple in this case be- 
cause only the exchange functional is significant whileas 
the correlation functional is small (In this paper, we shall 
ignore the correlation functional because it is absent for 
a single hydrogen atom.). We focus on the error of the 
work function predicted by the GGA or the X a method. 

When a is finite, the occupied single electron states are 
described by the Bloch state |k >= J d 3 ru k (r)|r > with 
\k x \, \k y \, \k z \ < ir /a. In the framework of Hartree Fock 
theory, u k = u^ F satisfy the Hartree Fock equation 



i_o „, , r „ .OXV), 



3 ,U™*(r>)u™(r)u™(r') 



r - r' 



= e HF (kKf (r), 



X (r) 



(1) 
(2) 



where 



(3) 



with Rj representing the lattice points. The Hartree Fock 
orbital energy can be also written as 

e HF (k) =< k|T + (7|k > + ]T(pp|kk) - (pk|kp) (4) 



with T = -V 2 /2 and 



(k 1 k 2 |k 3 k 4 ) = / d 3 r 

j J' - i- 

(5) 

Because the Bloch states |k > are extended, the Hartree 
Fock wavefunctions u p (r) will be not modified when one 
electron in |k > is removed. According to the Koop- 
man's theorem, the ionization potential associated with 
the state |k > is simply given by — enF(k). By adding 
the correction from the correlation energy, one obtains 
the exact ionization potential. Since the state at the 



X point k = (7r/a, 7r/a, tt/o) is at the top of the elec- 
tronic band, the work function is approximately given by 
-e HF (A). 

In the framework of GGA, u k = u GGA which is deter- 
mined by the Kohn-Sham equation 

(-iv 2 + [/(r)+ J d?v'n{v)/\v-v'\)u GGA (r) 
+ v xc {n{v))u GGA {v) = e GGA (k) U GGA (r) (6) 



with the exchange-correlation potential given by 

, , u SE xc [n] 
v xc (n(r)) = 



(7) 



5n(r) ' 

where E xc [n] is the exchange-correlation functional and 
n(r) = ]>>£GA> )M GGA (r) (8) 



is the total electron density. The GGA orbital energy 
can be also written as 

£GGA(k) = j d 3 ru GGA *(r)(-±V 2 )u GGA (r) 



+ 



Mr) + j 



, n(r') 



+ w xc (n(r))]n k (r) (9) 



with 



n k (r)=u GGA *(r)u GGA (r). (10) 
The total electronic energy of this system is given by 

E = ^ I d*vu GGA *{v)(- l -V*)u GGA {v) 

k J 

+ J d?vU{v)n{v) + X - J d 3 rd 3 r' n ^ M ^ ] 

+ E xc [n] + U NN , (11) 

where 

^ = E |B r ( 12 ) 

Next, we consider the N — 1 electrons system with the 
|k = X > state unoccupied. Since \X > is extended, 
the total electron density n'(r) for the N — 1 electrons 
system is approximately equal to n(r) in the large N 
limit. Thus the GGA solution u^ GA (r) is not modified 
and we self-consistently confirms that n' rts n. The total 
ground electron energy for the N — 1 system is given by 

E' = E/ ^^ GA *(r)(-iv 2 ) U GGA (r) 

+ /AC/(rK(r) + l|W^ 



+ E xc [n']+U NN , 



(13) 
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where means the summation over all k states except 
the X state and n'(r) = n(r) — njf(r). Since nx <C n in 
the large A*" limit, we obtain the work function 

W = E' — E = -e GGA (X) (14) 

which is just the orbital energy of the state \X >. Eq.(14) 
can be also obtained by the transition state method gen- 
erally used in the X a method. It is worthwhile to point 
out that for finite system or when the hole state is local- 
ized, Eq.(14) is no longer valid. 

From many calculations for atoms and small molec- 
ulars, we know that the orbital energy obtained in the 
GGA or the Xa method is significantly lower than the 
exact ionization potential. Therefore, one suspects that 
W is significantly underestimated by Eq.(14) in solids in 
which the valence electrons cannot be described by the 
plane wave states. In the following sections, this suspec- 
tion has been proved to be true. 



III. NON OVERLAPPING CASE 

In this section, we study the large a case in which the 
overlap between the atomic wavefunctions on different 
sites is negligible. In this case, one can immediately ob- 
tain the solution in the Hartree Fock theory and in the 
GGA or the X a method. Actually in the Hartree Fock 
theory, we have 

u k F 0) = 7= H e tk - R '0^ F (r - Rj) (15) 
* R, 

with <^ F (r) being the exact ground state wavefunction 
for a single hydrogen atom located at R = 0, and 

e HF (k) = -0.5. (16) 

Therefore, the work function predicted by the Hartree 
Fock approximation is W KF = 0.5. There is a correction 
to W from the correlation energy in the N — 1 electrons 
system corresponding to the hole polarization term. The 
magnitude of this correction is —37 /a 4 (We have ignore 
the Van der Waals interaction term which is very small.). 

On the other hand, in the GGA or the X Q method, we 
have 

^ GA (r) = ^^ e ^0 o GGA (r-R. 1 ) (17) 

with <f) GGA (r) being the ground state wavefunction in 
the GGA or the X a method for a single hydrogen atom 
located at R = and GGA (r) satisfying the following 
equation 

2 r J |r — r '| 

+ Vxc (n GGA (r)))0 GGA (r) = e GGA^ G GA (r); (lg) 



where 

n GGA (r) = \<f> GGA (r)f (19) 

with 

£ GGA (k) = e GGA = J d 3 r^ GGA *(r)(-iv 2 )0 GGA (r) 

+ [-l+Jd^^p^+v xc (n GGA (r))]n GGA (rl20) 

The above equation indicates that the work function W 
is indeed given by the orbital energy of the single atom 
in the GGA or the X Q method in this case. However, it 
is well known that the orbital energy of a single atom in 
the GGA or the X Q method is much less in magnitude 
than the corresponding ionization potential. For hydro- 
gen atom, e GGA is found to be —0.29 in X Q method, 
in which a is chosen such that the self interaction en- 
ergy is exactly canceled by the exchange functional E x [n] 
(Note that we have used the exact wavefunction <ft^ F (r) 
in Eq.(20) to calculate e GGA ). The difference between 
the work function given by the Hartree Fock theory and 
that given by the X Q method is (00|00)/3 « 0.208 with 

(00|00) = / ri 3 rrfV n ° F , (r)n ^ (r/) = 0.625. (21) 
J |r-r'| 

It should be pointed out that the difference between the 
exact e GGA evaluated by (/> GGA (r) and the above value 
is msig nificant@. In GGA, the magnitude of c GGA is 
even less than 0.29. Thus, the work function is remark- 
ably underestimated in this case by the GGA or the X Q 
method. 

We need to point out that the above conclusion is not 
only applicable to the hypothetic solid in which the spin 
polarized hydrogen atoms are located on a cubic lattice, 
but also applicable to the real crystalline solids, with spin 
up and spin down electrons, satisfying the following two 
conditions (i) The energy bands are either completely 
filled or completely empty, (ii) The lattice constant is suf- 
ficiently large such that the overlap of the wavefunctions 
between different lattice sites is negligible. This conclu- 
sion can be derived using the exactly same way as for our 
hypothetic solid. Many real solids for instance, the solids 
formed by Ne, Ar, Kr, Xe atoms at low temperatures and 
even ice (Note that there is polarization induced by the 
dipole moment by the nearby water molecular. But this 
effect is very small.) may satisfy the above two con- 
ditions. It should be noticed that for instance for the 
Ne solid, the work function is obtained by the Hartree 
Fock orbital energy together with the corrections from 
the electron correlation effect as well as the lattice relax- 
ation effect (i.e., the polaron effect). The magnitude of 
those corrections should be less than 0.1. However, the 
difference of the orbital energy between the HF theory 
and the GGA is estimated as 0.4. So the underestima- 
tion of the work function for the solid Ne given by the 
GGA is expected to be quite significant (about 0.3). 



4 



Although the ionization potential is remarkably un- 
derestimated in the non-overlapping case, the total elec- 
tronic energy E is still correctly given by the GGA or the 
X Q method when the band is completely filled. Since 



n GGA 



(r) 



,GGA 



(r - R t 



(22) 



and 



^(r-RO^V-R^-O (23) 
when Rj ^ Rj, we have 

£ = 7v|d 3 r0 o GGA *(r)(-iv 2 )0 o GGA (r) 
+ \r\ + \l ^^^ GGA (r) + E xc [n GGA ] 



-0.5JV. 



(24) 



Since the bandwidth is vanishing in the non- 
overlapping case, it may be interesting to consider the 
localized hole state. In this case, the ionization potential 
equals 0.5 — 37/a 4 and is correctly given by the GGA or 
the X Q method. Unfortunately, the GGA predicts that 
the energy of the extended hole state is lower than the 
energy of the localized hole state. 

It may be pedagogic to discuss the case that the band 
is only partially filled. Here, for simplicity, we shall only 
consider the infinite a case. Suppose that there are Z < 
N spin polarized electrons in the system. Assuming that 
the single particle states are all extended over the whole 
lattice. In the X Q theory, the total electronic energy is 
found to be 



E(Z,N) 



1 



-Z(l 



(00|00) z 
2 



(25) 



(Note that this equation is valid only when a — oo. When 
a is large but finite, there is an additional term C/a for 
E(Z, N) when Z ^ N.) In Eq.(2), we have used ^ F (^r) 
in the construction of Bloch states, where £ is a varia- 
tional scaling parameter which is used to minimize the 
total X a energy. Therefore, the total electronic energy 
is significantly overestimated when the band is only par- 
tially filled. Interestingly, when Z <C N, one may even 
obtain a lower value for E(Z,N) by assuming all single 
particle states are extended only over five lattice sites. 
According to Eq.(25), one may obtain the ionization po- 
tential for Z electrons 

I(Z, N) = E(Z - 1,N) — E(Z, N) 



-5" 

[1- 



(00|00) Z 
N 
Z 



2 

(00|00) 



[3 



N 



-(-) 1/3 ]] 
\ N J JJ 

3 y N J i 



(26) 



Notice that when the band is half filled Z = N/2, we 
have / = 0.516 which is quite close to the exact value 0.5 



It may be necessary to mention that this deficiency in 
the GGA or the X Q method is also found [l(| HH, El m 
some molecular systems at the dissociation limit, such 
as Hj, Fj, NaCl, et.al. Essentially, both failures of the 
GGA or the X Q method for the solids at large a and 
for moleculars at the dissociation limit are caused by the 
same reason. The authors in Ref. [HI EH, E2| attributed 
this type of failure to the so called self-interaction error. 
However, we deem that this viewpoint is inappropriate 
because the self-interaction error is zero in our system, 
yet the magnitude of the total electronic energy is still 
significantly overestimated by the X a method. This fact 
also suggests that the self-interaction-corrected local den- 
sity approximation method (l3l . [3] is not helpful to im- 
prove the result. 

We notice as well that the overestimation of the magni- 
tude of the electronic energy by the X Q method is more 
pronounced in our system than in diatomic molecular 
cases, because the electronic states are much more ex- 
tended in our case. For instance, consider the ioniza- 
tion potential of two spin polarized hydrogen atoms at 
the fixed separation distance a. When a is large and 
the atomic wavefunction overlap between two hydrogen 
atoms is negligible small, we can obtain that the error of 
the ionization potential given by the X a method 



SI 



(00|00)[(|) 4 / 3 



1 



0.092 



1 

4^' 



(27) 



where the polarization term has been neglected because 
it is proportional to 1/a 4 . This error is significantly lower 
than that in solid case. For instance, when a = 5, this 
number is only 0.04 which is much less than 0.21 in the 
solid case. This fact explains why the GGA or the X Q 
method are still successful for small molecular near equi- 
librium geometry. 



IV. FINITE OVERLAP CASE 

In the last section, we have discussed the large a case in 
which the overlap of the wavefunctions on different sites 
is negligible. However, the GGA or the X a method are 
usually applied to real solids of finite bandwidth. When 
a is sufficiently small, the electrons in our system can be 
described by the plane wave. Therefore, one may specu- 
late that the GGA becomes valid for not large a case in 
our system. On the other hand, for large a, the atomic 
nature still remains and it can be expected that the er- 
ror of the work function predicted by the GGA or the 
X Q method is still significant even when the bandwidth 
is finite. So it is highly desirable to see how this error 
decreases when a decreases. 

When a is finite and the bandwidth is not zero, it is not 
easy to calculate the electronic energy accurately even 
within the framework of the HF theory or the GGA. In- 
stead, we shall take several approximations which can 
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x k GA ( r ) an d calcu- 



greatly simplify the calculations. 

At first, we assume uj^ F (r) i 
late the difference of the work function predicted by the 
Hartree Fock theory and the GGA or the X Q method. 
When the hole state is extended, straightforward calcu- 
lation yields 



AW = Whf - W GGA = Ai - A 2 (28) 



with 



A 2 



Ax = £(kX|Xk), 



d 3 ru xc (n(r))n x (r) 



(29) 



(30) 



where n and nx are the total electron density and the 
electron density of the state \X >, respectively, Whf and 
Wqga are the work function obtained by the Hartree 
Fock theory and the GGA, respectively. Note that this 
assumption has been tested in atomic calculations. The 
error from this assumption is found to be insignificant [|| . 
Secondly, we assume 



"k(r) 



:HF 

'0 



(r - R, 



with the normalization constant 

f 



C(k) = -= 
where the overlap is defined as 



(31) 



(32) 



(33) 



This assumption is valid in the tight binding limit. While 
for small a and large bandwidth, different |k > states are 
made of different <f> and this assumption becomes invalid. 
However, since both the quantities < kA|Ak > and v xc 
do not depend on Uk sensitively, we expect that this as- 
sumption is valid at least when a is not small. Under 
this assumption, the work function difference AW can 
be evaluated numerically. 

Thirdly, we shall use the X a exchange functional with 
the parameter a — «h selected to cancel the self- 
interaction exactly in the atomic limit. So we have 

nr. 1/3 f 

E xc [n(r)} = J dW/3 (r) . (34) 

Since an is greater than the corresponding value in the 
electron gas model, and furthermore our system will re- 
duce to the electron gas model in the small a limit, wc 
expect that this exchange functional will slightly overes- 
timate the exchange energy in the finite a case. Table. 1 
lists the numerical value of Ai, A2 and AW for vari- 
ous values of a. From this table, one can see that when 



a 


3.0 


4.0 


5.0 


6.0 


7.0 


00 


Ax 


0.859 


0.738 


0.685 


0.656 


0.639 


0.625 


A 2 


0.716 


0.577 


0.504 


0.465 


0.442 


0.417 


AW 


0.143 


0.161 


0.181 


0.191 


0.197 


0.208 



TABLE I: The calculated results for Ai,A 2 ,AW at various 
values of a. 



a > 5, AW significantly deviates from zero (Note that 
when a = 3 and a — A, we expect that the second approx- 
imation becomes invalid. For instance when a = 4, the 
bandwidth is estimated as 0.5 and the electronic states 
cannot be described by Eq.(31).). Even after subtract- 
ing the correction from the correlation energy contribu- 
tion 37/a 4 , the underestimation of the work function by 
this X Q method is still about 0.12 at a = 5. On the 
other hand, the band width is about 0.25 when a = 5. 
Therefore, we can conclude that the work function will 
be underestimated by the X Q method when (i) the va- 
lence bands are completely filled and (ii) the bandwidth 
of the valence bands is not large. 

We have also used the generalized gradient approxi- 
mation [i| exchange functional to calculate A2. In GGA, 
the exchange functional is expressed as [4| 



E xc [n] = Ajd 
with A = -(3/4)(6/tt) 1 /3 ; 



3 rn 4/3 



F(t) 



(35) 



F(t) = (l + 0.021326i + 0.0037909i 2 + 0.000000891i 3 ) 1/15 

(36) 

and 



t = 

And A2 is evaluated by 
A 2 = A 



(Vn) 2 

77,8/3 



(37) 



-F'(t)[ 



2V77-V77X 8 (Vt7,) 2 ti. Xi 



,4/3 



,7/3 



(38) 



As expected, A2 evaluated by this formula is indeed 
slightly less than that by the X a method. For instance, 
when a — 4, A2 is found to be 0.548 in GGA which is 
slightly less than the value 0.577 given by X a method. 



V. DISCUSSION 

In this paper, we have studied the hypothetic system 
composed of the hydrogen atoms with spin polarized elec- 
trons on a cubic lattice. We found that for the case of 
narrow bandwidth and completely filled band, the work 
function is significantly underestimated by the GGA or 
the X a method. This is attributed to the significant lower 
energy of the extended hole state obtained by the GGA or 
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the X Q method. Note that this problem cannot be cured 
by assuming that the hole state is localized. Actually, 
the GGA or the X Q method predicts a significant lower 
energy for the extended hole state than for the localized 
hole state even in the zero bandwidth limit. 

Although we have only calculated a hypothetic system, 
our results may lead us to speculate the validness of the 
GGA or the X a method in the calculations of the work 
function of the insulators with narrow valence bands in 
many ionic crystals, as well as Ne,Ar,Kr,Xe solids (as- 
suming the lattice constant fixed.). Furthermore, it is 
not difficult to see that the GGA or the X a method be- 
comes invalid in calculating the ionization potential for 
the extended, including noncrystalline, systems, where 
the valence electrons cannot be described by the plane 
waves. As we have already observed, the GGA or the 
X Q method would predict that the hole state is extended 
over many different regions which are distant from each 
other even in inhomogeneous systems. Such prediction 
differs from the commonsense. 



In section III, we have also pointed out that the mag- 
nitude of the total electronic energy is significantly over- 
estimated by the GGA in the large a limit when the 
band is not completely filled. Actually, Zhang and Yang 
and other groups [l(J [HI G3 also found that H^Fg , 
NaCl et.al. would not be dissociated correctly in the 
GGA in the large distance limit. This deficiency is more 
pronounced in solids because the electronic states are 
more extended. It should be emphasized that the self- 
interaction of the electronic states in our system is negli- 
gibly small. Therefore, this deficiency cannot be removed 
by the so called self-interaction-corrected local density 
approximation method(l3l. fl4|. 

In summary, we speculate that the exact exchange- 
correlation functional is highly nonlocal when the elec- 
tron density n(r) varies in space and the system is ex- 
tended in the density functional theory When the 
density of the valence electrons varies greatly in space 
and the system is extended, the GGA and the X Q method 
may be unsuitable to describe the system properly. 
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